function [f, df, d2f] = functest1(x, alpha)
    % Calcola la funzione nel punto x
    f = x.^2 - alpha;
    
    if nargout > 1
        % Derivata prima
        df = 2 * x;
        if nargout > 2
            % Derivata seconda
            d2f = 2 * ones(size(x));
        end
    end
end
